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ABSTRACT 


The effects of fuel burnup, fission product poisoning, and hydrogen 
moderator density variation on reactivity in the AGN-201 reactor are 
considered. A modified one-group perturbation theory is developed and 
applied to changes in parameters resulting from a change in hydrogen 
moderator density. An equation for the feacwiveey change is obtained 
for three models: a bare cylindrical core, a bare core using extrapo- 
lated dimensions, and a reflected core. These three equations are then 
used to predict values of the reactivity increase resulting from inter- 
changing a new 3/4 inch thick fuel disk with comparable fuel disks 
presently in the core. 

rhe results obtained by a digital computer solution of the re- i” 
activity equations reveal that the increase in eactiuley varies from 
0.4392 to 0.707%, depending upon the core model daa position of the old 
disk within the core, Because the license of the: Havel Postgraduate 
School deus not permit a aide of excess reactivity above 0.40%, it is Fe 
concluded that a simple interchange of disks in this manner would pro- 


duce too large a value of excess reactivity. 
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ENGLISH 


NOMENCLATURE 


SYMBOLS 


molecular weight of polyethylene 
atomic weight of uranium 

buckling of the reactor system 
extrapolation distance 

diffusion coefficient of core 
diffusion coefficient of reflector 
thermal utilization 

non-1/v factor 


distance of nearest edge of fuel disk above or below the zero 
ordinate of the z-axis of the cylindrical fuel core 


height of fuel core 

extrapolated height of fuel core 

effective resonance integral 

multiplication factor for an infinite reactor system 
diffusion length in core 

diffusion length in eer tacter 

perturbation operator 

adjoint of M 

perturbation operator of perturbed reactor system 
mass of fuel per disk 

mass of polyethylene 

atom density 

Avogadro's number 


atom density at time equal to zero 





P 


resonance escape probability 


Pm = promethium concentration 


PWR = reactor power 


R 


R! 


radius of fuel core 


extrapolated radius of fuel core 


samarium concentration 


equilibrium samarium concentration 


time 


reflector thickness 


effective neutron temperature 


moderator temperature 


radial reflector thickness 


end reflector thickness 


disk volume 


centerline distance from a disk to the zero ordinate of the 
z-axis of the fuel core 


GREEK SYMBOLS 


bp 
3 
sh 
Sr 
32 


An 


Ke 


lt 


fission-product poison yield of promethium 
reflector savings 

disk thickness 

radial reflector savings 

end reflector savings 


fraction of total change in core moderator mass occurring 
in a fuel disk 


fast fission factor 
net neutrons available after fission 


decay constant of promethium 


Ne 


§ a3s 


a MW 
° a 


Pp Be 


7 


~*~ BR 


transport mean free path 
average number of neutrons released per fission 
average logarithmic energy decrement per collision 


reactivity 


density of polyethylene 


density of uranium 

density of U-235 

density of U-238 

microscopic absorption cross-section 
microscopic fission cross-section 
microscopic scattering cross-section 
macroscopic absorption. cross-section 
macroscopic fission cross-section 
macroscopic scattering cross-section 
Fermi Age | | 

neutron flux 

average neutron flux 

peak neutron flux 

neutron flux of perturbed reactor system 


importance function 
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I. INTRODUCTION 


The AGN-201 nuclear reactor is a small, safe, low-power, self- 
contained reactor designed and produced by Aerojet+General Nuc leonics 
Corporation of San Ramon, California. It is designed primarily for use 
in education and research not requiring a high neutron flux, but featur- 
ing maximum safety, high analytical sensitivity, and ready accessibility 
of components. The fuel core of the reactor is composed of ten circular 
fuel disks of varying thicknesses, piled on top of one another to form 
a cylindrical shaped core. The fuel disks are made of polyethylene and 
have uranium fuel homogeneously interspersed in them. Appendix A con-. 
tains a general physical Sesneieeien of the overall reactor and also a 
detailed discussion of its fuel core. © | 

The AGN-201 nuclear reactor located at the Naval Postgraduate School, 
Monterey, California, has been. in licensed raat Gn since April 29, 1957. | 
From that date, until January 31, 1962, the reactor was operated at a 
maximum power level of 100 milliwatts. During that time period no de- 
crease in axceue reactivity was noted. Thereafter, the license of the 
Naval Postgraduate School was amended to permit continuous operation up 
to 20 watts and intermittent operation at power lavels io to 1000 watts. 
On February 3, 1962, the excess reactivity of the reactor was measured 
and recorded as 0.354%. From that time until February 26, 1969, the 
value of the excess reactivity has steadily decreased to a value of 0.281% 
or a drop of 0.073%. 

| Several reasons which would account for the observed loss in re- 
activity have been postulated. The first cause of this loss could be 


due to fuel burnup. Since the reactor is highly enriched with uranium-235, 





(approximately 20%), the decrease in reactivity would be proportional to 
the fraction of U-235 depleted by fission. A second cause of reactivity 
drop could be the buildup of reactor poisons due to fission of the fuel. 
These fission fede and their decay products have large neutron eae 
tion cross-sections, If they are preduesa in appreciable amounts they 
could effectively alter the neutron balance in the reactor by affecting 
neutron multiplication. Thirdly, the loss in reactivity of the reactor 
might be attributed to a loss in apdeveude within the saatae core. This 
problem is unique to this type of aactok, Polyethylene, which is the 
neutron moderator, when subjected to radiation, tends to decompose. This 
decomposition is accompanied by the liberation of hydrogen gas. The ef- 
fect of this phenomenon on the reactivity is especially pronounced be- 
cause the hydrogen is bled off from the fuel core at varying time inter- 
vals which reduces the density of the moderator within the core. 

In September, 1968, the Naval Postgraduate School received a new 
fuel core disk for the purpose of adding it to the core in some manner 
to attempt to raise the value of the excess reactivity of the reactor. 
However, the limitation imposed upon this addition is that the resulting 
value of the excess reactivity of the reactor must not exceed 0.40%. The 
new disk has a known fuel enrichment and a higher polyethylene density 
than any of the disks presently in — core, The properties of the new 
disk are presented in Appendix A. . 

The objective of this thesis is to determine the change in reactivity 
which would result from interchanging the new fuel disk with a disk of 
comparable size currently in the reactor. To accomplish this objective, 
first an analysis is made of the three possible causes of the reactivity 


loss due to reactor operation to determine which have significant effects, 
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and what changes in reactor properties result from these effects. Second, 
using a modified one-group perturbation theory, an equation for the re- 
pbigety as a function of these changes is derived for three models des- 
cribing the reactor core. These models treat the core separately as a 
bare cylinder, a bare cylinder having extrapolated Soueneienn, and a re- 
flected cylinder. Third, a digital computer program is used to solve the 
reactivity equations using assumed changes in reactor pRouere tee sane 
the actual value of reactivity drop due to operation is known, the core 
properties at present can be found. Using these results and ye proper- 
ties of the new fuel disk, the reactivity equations derived from perturba- 
tion theory are solved by a second computer program to yield values of the 


excess reactivity caused by the perturbation of replacing an old disk with 


_the new disk. 











II. VARIABLES WHICH AFFECT REACTIVITY 


im FUEL BURNUP 

The first of the three probable causes of decrease in excess reactiv- 
“tty toe considered is that of fuel burnup. From published literature 
weoneeratne the AGN=-201 reactor fy, it is predicted that continuous re- 
aneee operation (24 hours a day, 7 days per week) at a 100 milliwatt power 
“Pavel woatd result only ina burnup rate of approximately 30 micrograms/ 
. cage. On the basis of this information it can probably be anticipated 
that even ae donee power levels, the contribution of fuel burnup to a 
upeactivits drop will be slight, if any. 

Since. the reactor is fueled with only a single fissile substance, 
U-235, and there are no fertile materials present, (it is assumed that 
due to the high enrichment. of U-235 there is no compensation for fuel 
depletion through the production of fissile Pu-239 from neutron absorp- 
tion by U-238), the depletion rate of the fuel is governed by the equa- 
tion 


aN S 
ae = -NDS (1) 


where N is the atom density of the fuel, F, is the thermal absorp- 
tion cross-section, and @ is the neutron flux. Rearranging and 


integrating (1) yields 


N [Eat at 
= @ 


No 


t ~~ — = 
The integral exponent, f G, %t dt , can be expressed as GGT , 
& 


— 


where @ is an average value of the flux over time t. Then 


 Nombers in brackets refer to Bibliography. 
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N. geet ie 
From recorded data contained in the operator's log kept at the reactor 
facility, as of December 9, 1968, the reactor had been operated approxi- 
mately 8.69 x 10° watt-min. Although operation occurred at seane levels 
varying from 0.1 siete-ne 1000 watts, 20 watts was most frequent, and 
thepetoes Was assumed £6:'6e a representative average value. “This gives 


an average time of operation of 


t = 2.61 x 10° sec. 
An equation given by Lamarsh [5] relating reactor power and neutron flux, 


assuming a recoverable energy of 200 Mev per U-235. fission is 


=. CPWRITe) x IOS 
PE CIAK MENGE) @) 


where PWR is the reactor power, T is the effective neutron tempera~ 


eff 


ture, M, is the mass of fuel, and Be is the non- Yy factor. Since 
PWR = 20 x 107° megawatts 
Tore = 340° K- 
g.(T) = 0.9269 

and My = 0:667 Ke 


evaluation of equation (3) yields 


= 8,3 x 10° Saukione [en bes 


Sw 
i 


and the quantity 
6 


Bt = 1.16 x 10° (4) 
where 
®> saan VES T) 
| | a Tae VT Ta 9.0 
Substituting (4) into (2), 
N = Chie AIO" *) 
Ne » @ = 0.9999988 
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or No 1. 


Therefore it is concluded that fuel burnup is negligible at the assumed 


average power level of 20 watts. . , a 


B. FISSION PRODUCT POISONING 
As aravionts mentioned, reactor poisons, having large neutron absorp- 
tion ehogesacctions, have an adverse effect on reactivity. Their chief 
Mneldenee on reactivity is through the multiplication factor by decreas- 
ing the thermal utilization. In discussing reactor poisons in the AGN- 
201, it is felt that in order to affect reactivity, a poison produced in 
| this particular reactor must meet three criteria, (1) It must have a 
high cross-section. (2) It must havea solid form. (3) It must be “ 
stable. Although as many as 200 species are apt to be produced during Cs 
the fission process, only two have aes enough cross-sections and are a 2 
produced in sufficient quantity to warrant attention; namely Xe-135 and 
Sm-149. 
In aaNet these two, it is seen that Xe-135 meets neither condi- 
tion (2) nor (3) above. Being an inert gas, it would be removed from 
the reactor at the time when the hydrogen buildup is bled off from the 
core. Thus the decrease in reactivity would be a function of time, being. 
zero at a time equal to zero, buildup to some value, and then return to 
zero at the time of the hydrogen and poison removal. However, this is 
contrary to observed evidence of reactivity behavior. In addition, Xe- 


135 is not stable but is an intermediate isotope in the chain 
198 195 136 135 135 
Se a | Ne «eG ae (86 (stable). 


Because it is not a stable isotope, any xenon poison produced during i 


operation would eventually decay after reactor shutdown, resulting in 








no long term effect on the reactivity. Hence, Xe-135 can be effectively 
ruled out as a cause of long term reactivity drop. 


$m-149, on the other hand, meets all three requirements. It has a 





thermal absorption cross-section equal to 58,500 barns at 20°C. It is a 


eee tae rare earth metal and is the stable decay product of the chain 


m9 44 
Na — Pen —- Sm fo (stable). 


le , The neodynium decays relatively quickly compared to the promethium so that 


for practical purposes the Pm-149 may be assumed to be produced directly 





as a result of fission. The concentration of Pm-149 is then governed by 


“the equation 


APm 
at 





where Pm is the promethium concentration, So is the fission-product 


a 


poison yield for promethiun, ae is the fission cross~section of the 


« 


fuel, @ is the neutron flux; and Ap is the decay constant for prome- 
dium: 

The rate of change of the concentration of Sm-149 can be waptenued 
as 


Ge = AoPm - ag S 6) 


since samarium is stable and can only be lost through neutron absorption. 


In equation (6) S is the samarium concentration and C is the absorp- 


ces excep aber yA aor cintr ee ie prawn rfp tbeinata inate gt oat stom Daay serio mark 


tion cross-section of samarium. Combining (5) and (6) and integrating 


yields an expression for Sm-149 concentration as a function of time. 


& 


, ~ Aet -os gt : 
. s+ TEE (l-e Xie oP") (7) 
I 
| 
| i5 











tt 


| Since Yep = 0.0113 


S_ = 0.0646 cm 
5.85 x 410°°° ene . 


3,56 x 107 ade 


~ 9 
D> Lan 
tt a 


and from the previous section it was assumed that 


2.61 x 10° sec 


t 


% 





8.3 x 10° neutrons/cm”-sec, 
Then 
1.313 x 10" eur? 


S 
It can be seen from equation (7) that as t —*® , S would no longer be 
a function of time but would approach a constant "equilibrium concentra- 


tion”. For purposes of comparison let t —*© in (7). Therefore, 


Ye Ze 16 


Ses - “Rr = 1.178 x 10 ei "s : 





“Thus it ia sean thet the concentration at present is less ‘San ene: 

equilib eta value by a factor of approximately 10", The macroscopic 
cross-section of the samarium is | 

a Wen te 8 = Sa = 7.68 x 108 em. 


The change in reactivity due to samarium poisoning can be determined 


from the equation [5] 


§ = WAY 
wEP 
where 
DW = 2.43 
and . ps € m= il. 
Thus, 
§ = -0.48 x 107". i 


Using the same equations with the equilibrium concentration of Sm-149 
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In fact, the molecular weight of commercial polyethylene is between 


ieits 
§ = -0.439% 

As can be seen, there is a difference of a factor of 10" between the 
reactivity values. Further, using the same equations to solve for the 
time required to reach equilibrium concentration yields a value of 8.11 
x 10’ sec. This is a factor of 40 times larger than the time of re-— 
actor operation to present. 

From this analysis it can be concluded that samarium poisoning does 
exist in the reactor, but it has not reached its equilibrium value. 
Due to the relatively short period of operating time and the low flux level 


produced during operation of the reactor, the decrease in reactivity due 


_to samarium is several orders of magnitude smaller than the observed value 


of 0.073%. Consequently the effect of samarium will be neglected.  , 


C. MODERATOR CONCENTRATION | 

The third possible seigecee reactivity drop is that of loss of modeta- 
tor. As already noted, the moderating material in the AGN-201 reactor is 
polyethylene. Chemically, the aatacalag formula of polyethylene is (CHa 
indicating that tt de composed of long chain (high polymer) molecules 
formed by many (CH) groups. | 

| ~CH, CH, CH, CH, ~ 

18,000 and 22,000. From a nuclear standpoint, polyethylene is a long 
séries of carbon atoms connected in chains with bare protons (hydrogen 
nuclei) attached to them. When high energy wenbeehe nes introduced into 


a polyethylene medium, it is possible for.a neutron to collide with the 





proton and break the proton away from the main chain. When this happens, 


cross-linking occurs, whereby new bonds are formed between adjacent high 








polymer molecules accompanied by the evolution of a hydrogen molecule. 


' CH, CH, CH, CH, = _ CH 


CH, CH, CH, CH,.CH. CH, = 


2°°2 2 hes 
ANd | HL, 2 
~CH, CH, CH, CH. - . ORES CHS CH SHG S. * 
This is referred to as unsaturation of the polyethylene. The main point | 
to be noted is that when polyethylene is irradiated, gas is evolved, the 
majority of which is hydrogen, along with slight traces of methanes. The 
assumption made, upon all that follows in this thesis ie thet due to 
hydrogen evolution and subsequent removal of the hydrogen frou: tha ‘fuel 
core, the mass of the polyethylene moderator is decreased and that this 
decrease is equal in magnitude to the mass of hydrogen removed. 
‘From research cotidiicted with the AGN-201 in March ,1962, shortly after 
the teuctor license was aneuded to permit higher power operation, it was 
noted that gas evolution from the core occurs primarily at power levels 
., between 200 and 1000 watts. [7] This gas evolution takes place.slowly | ; 
and aga a period of about two days following operation. No gas evolu- | 


tion was observed in relatively short periods of operation at 20 watts. 


A sampling of the gas revealed traces of an inert gas mixed with the 





nydvoses: the decay rate of which corresponded to xenon. The fact that 
the beginning of the observed decrease in excess reactivity coincided with 
the start of the reactor's operating at high power levels tends to con- 
firm that hydrogen evolution is the cause of the reactivity drop. 
| In the theory that follows, the decrease in density of the poly- 
ethylene moderator manifests itself in two ways. | First, it. appears . 
_ directly in the moderator density term, S a Secondly, it causes a change : 
in the fission and absorption cross-sections, bs and T es of materials 
in the reactor. Because the amount of moderator in the fuel core is 


decreased, the neutrons are not as well moderated and consequently the 
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average energy of the neutrons is at a higher energy value. When this oc- 
curs, the cross-sections of the materials in the reactor, which are 7 
energy dependent, are no longer those corresponding to 0.025 ev, but rather 
those which correspond to a slightly higher energy. 

The energy at the peak of the energy curve is related to a tempera- 
ture known as the effective neutron temperature. The neutron temperature 


is related to the density of the polyethylene in the following manner. 


- Coveyou [9] has established the relation 


| a ae Ty (1 +.0.46 A) 
ons | A = Fz a? 2.) 
and | gin S?No 
m rp 


for moderators with mass from 1-to 25 and Q from 0 to 1.0. For heavy 


moderators, Cohen [9] gives 


Tore = Th (1 + 0.34 3 
From a survey of research done in the field of polyethylene, Cooke [3] 
concluded that the effective neutron temperature, Tore? of. neutrons in 
polyethylene was between 336 and 347° K. Since it is known for the re- 


actor that A = 0.2931 and Ts = 293° K, and To¢¢ Wa8 assumed to be 340° K, 


ff 
the numerical coefficient in the expression for Tope Was modified to a 
value of 0.5459 to fit the polyethylene data. The resulting equation is 
thus, 


Tope ~ Ty (1 + 0.5459 G4, (8) 


From the equation linking the density of polyethylene to the effective 
neutron temperature, the neutron temperature can in turn ‘be used to find 


the cross-section at that temperature through the following equation, 
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Thus, from these equations it can be seen that the effect of decreasing 
‘the polyethylene density is to increase the neutron temperature which 
decreases the neutron cross-section of materials in the reactor.- 

From work done with the reactor in June, 1963, based on test runs 
of 1000 watt-min, the rate of hydrogen evolution from the fuel core. was 
found to be about 0.72 moles/KW-hr. [8] The number of kilowatt-hours 
of operation from February, 1962,to February, 1969, the period during 
which the drop in excess ere has occurred, is known to be about 
16.07 KW-hr. Therefore, the number of moles, and consequently the weight 
of hydrogen gueivea oan be approximated. . This results in a value of — 
hydrogen weight of 23.14 gm which when equated to a chanverts polyethy- 
lene mass in the core gives a value which can be used to check against 


that predicted by perturbation theory. 


D. ° ‘CONCLUDING STATEMENTS 

- The conclusions drawn from the preceding sections can be summarized 
as follows. Fuel burnup can be neglected as a cause of reactivity de- 
crease on the grounds that the time of operation of the reactor has not 
been long enough nor is the assumed average level of operation at 20 
watts high enough to cause any significant depletion of the U-+235 fuel 


; in the reactor core. The two reactor poisons most likely to be produced 





during operation of the AGN-201, xenon and samarium, cannot account for 
‘the drop in excess reactivity. The Xe-135 is either withdrawn from. the 


reactor simultaneously with the hydrogen during bleedoff or decays to a 


more stable isotope during the frequent periods of reactor shutdown. The 
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extremely small concentration of Sm-149 present in the reactor does pro- 


duce a decrease in the excess reactivity but its magnitude is several 


orders sna ller than the observed decrease and is negligible. The de- 


crease in polyethylene wcieeaces due to hydrogen evolution resulting from 


‘neutron irradiation cannot be ruled out as a cause of reactivity drop. 


It ig assumed that this cause is solely responsible for the decrease in 
reactivity. Also, observed evidence seems to indicate that this is the 
most probable cause of the reactivity decrease. The change in polyethylene 


denaity will appear in the theoretical equations that follow both as a 


direct density change and also as changes in neutron cross-sections of 


reactor datevtale: This is because the cross-sections are functions ‘of 


neutron temperature which in turn is a function of moderator density. 





ITI. PERTURBATION THEORY 


A. INTRODUCTION 

To analyze the effect on the behavior of a nuclear reactor caused by 
‘small non-uniform changes such as those discussed anti preceding pata- 
graphs, a method known ag Aatkunpation eH eoey: has been daveione’: If the 
problem, which tias been sabiiaed’ involving the AGN-201 were of such a 
nature that uniform aienaeaiacourted throughout the reactor, then kha re- 
activity could simply be found by recomputing new reactor constants and 
recalculating a oe Eatotaeation factor to find the resulting reactivity. 
When non-uniform changes thie the reactor occur, as is usually the case 
in practice, then perturbation theory must be employed. Examples eee 
uniform perturbations include insertion or withdrawl of control rods into 
the reactor, non-uniform fied burnup, aoa ml Poea fvesion poison accumu= 
lation or as it is believed in this case, sioneunttone removal of modera- 
tor from the core. The assumption upon which perturbation theory is based 


is that the perturbations or changes within the reactor are small and that 





the flux is not distorted substantially in the region of the perturba~— 
tion. Since this is assumed to be the case in analyzing the AGN-201, 
perturbation theory provides a convenient method of handling an other- 


wise difficult problem. 


B. . MODIFIED CRITICALITY EQUATION 
Perturbation theory is referred to as one-group, two-group, or many- 
group depending upon the type of theory used to describe the energy of the 


neutrons in a critical reactor. The simplest is one-group and will be used ix 


here. The criticality equation describing a one-group reactor is 








where Ke = ut 


Strictly speaking, this equation is used in a one-group perturbation 
analysis. However, it was felt that fappoved accuracy of the éauite. 
cold be obtained by using the Fermi Age-Diffusion model to include neu- 
trons in the slowing down region. The usual fg of the criticality 


equation in this case is 


= B*A, 
Ko 6 
Teter * | - 
- 08, 
where Ka* 7 pf ; € is the fast non-leakage probability, 


and Creute*y” is the slow non-leakage probability term. However, the 
method used to calculate a, (the Fermi Age equation), is not valid in 
the slowing deais region for the eee se an hydrogeneous moderator be-= 
cause the assumption of many collisions by the neutrons does not hold. 
If the moderator contains hydrogen to any degree, it ti: possibile for 

the neutrons to become thermalized in a few or perhaps even a single 
collision. Weinberg and Wigner [9], have developed a single collision ~ 
kernel to account for this type of moderation. The result is that the 


usual fast non-leakage term is replaced by 
Ze ~| g 
= tan”! HE) (0) 


where Zs is the cross-section for this type of interaction and is 
ae. ] 
iven as 2.7 —— 
: OY RSG.) Ye 


where Cs. is the experimental value of the Fermi Age. Substituting 


(10) into equation (9) gives 


Keo (3) tan7'C#) a 


pr bes 


Bs 
fd 











This is the form of the modified one-group criticality equation to. which 


strict one-group perturbation theory will be applied. 


c. DEVELOPMENT OF PERTURBATION THEORY 
In this section, one ~group perturbation theory is employed to derive % 
. an expression for a change in reactivity from changes in various nuclear 
parameters. ‘ 
Let Ce = ten” -z) 
The modified one-group criticality equation for an hydrogenous reactor 
deve Loped above becomes | 
1+ L8G? 


- This equation is the direct result of a modified diffusion equation of | 


= | 


the following form. 

VEPCE, P= LH +divDgrad PO — (11) 
where VEPCI.S _ is the source term, Za is the absorption term 
and div D grad @ is the leakage term. Equation (11) can be written as 


MZ=z0 
where M is the operator . 
| M= VEPCE,- Zatdw D grad 


Referring to the discussion of the adjoint operator in Appendix B, the 





adjoint equation, which will be of use later in the development but shiek 
is presented here, is 
Mt¥=0 | 

M+ is the adjoint of M and the function ¥ is referred to as the im- 
portance function. For one-group Saeeathabion cheony, M is self-adjoint 
and & and * are everywhere proportional. Terms in sauak ten (11) which 


are variables are P, C, Ze s Za , and D. Considering first perturbations 


24 ! 





or changes in P, C, 2, and Sa > 


pe 


ll 


P+ §P 
c+ $c 

Ze = Z,+ 82; 
Sa’ = Bot lq 


The unprimed terms refer to unperturbed quantities, the delta terms are 


c? 


magnitudes of the changes or perturbations, and ‘the primed terms. refer 
to the resulting perturbed quantities. To return the reactor to critical, 
a change must be made in WY . This is not possible in vegitey since YU 
is a physical constant. However, the concept of compensation in this man- 
ner ig useful in the analysis, so that 

v's vt Sv 
For the perturbed reactor 

M'g’=0 
where 7 : 
M! = (+ SDKEXKP + SPLC+ SCX Let $Z,) - C Lat Fla) raiv D grad 
Neglecting higher order terms, since for perturbation ‘theory to be. valid 
“ne perturbations must be small, we get | 
M' = DEPCE, * SVEPCZ, + SPVECE, + IC UEP TZ, 

+ $2, VEPC - Zy* Bdq + div D grad (12) 
Substituting the quantity M into (12) gives 
M's M+ svercd, + SPECS, + 82, VEPC - TEq 

This may also be written as 


Mie M+ Ag 





where 


Ag? SVEPCE, + SPUECE, +32,VE PC - B24 








Therefore C(M+A,) 0’? ° a a (13) 
Mtr eo (14) 


Multiplying. (13) by Y , (14) by @ , subtracting (14) from (13) and 
integrating the result over the volume of the fuel core yields 


of vO + Ae" dv - fy BM*pdveo 


-Expanding the first integral, 


~ feme'av+ Sf vag ave “Sy oH av 6 


Combining the first and third integrals, 
wt 7? ? + ae z 
St, Ceme’-@'mtvyav eS, #A.@'aV2 0 (a5) 
However, from Appendix B, by definition of the adjoint operator, 


JS, CeMa’= a’ Mtv) av = 0 


Therefore, (15) becomes 


J A,B aVEO 
Substituting the expression for Ag gives 
Sy PCSVEPRCE, + SPVECE, +H SCVEP Zp + 32g VEPC - 3Eq) @' dV = O 
Solving for y in since J = - yy $84 
Jy YCEPECEe + SC EP + BZ¢EPC - “TT 
SL, # CE PC Ze) w’ AV 


Since it is assumed that the perturbation is small, 
o =D 

and because the one-group operator is self-adjoint, 
peg | 


and the expression for J becomes 


_ [ucarecde * Waray + 38+€°C)- 32.] ¢*av 


(16) 
ff (Cwercay) o* AV 
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Next, consider a ones in D daly, 


D'= D + 2D 
Then 
7 : 
Mos (u+tvXePCEs) + aiv CD+SD grad 
Following the analysis. previously outlined, the operator for the per=.. 
. turbed system can be written as 


M's M+A, 


where the perturbation operator is 


A, * SvEercd, + div 30 grad 


This results in a reactivity of 


fe Lsdi s0eraav" 
ty CvEPCS,)@* av 
By employing vector identities, 


Jy #aiv Waraagav = Sf aiv #sDgrad@)av~ -xocow yay 


However, since the flux vanishes at the reactor core euraece 


(17) 


Sy div (esd grad Dav = Sy BO grad wei ahs 
Hence, 7 
J, Bdiv 8D graag ave ~S) sp cv)? av 
Therefore, (17) can be expressed as . 
. Jy SDCV8)* AV 
Sy (Ve PCE,) Sav 


Changes in P, C, a, Zo and D occurring simultaneously produce a re- 


(18) 


activity represented by the sum of (16) and (18) or 


e: J flvcrecss + 3CEPS, + $2p€0C)-3%q]8"- DVB) FAV. 


S, (VEPCE,) Bray 
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D. APPLICATION OF PERTURBATION THEORY 
1. Bare Cylindrical Geometry 

In this section the expression derived for the change in re- 
activity is solved for the: case in which the fuel core is unreflected . 
dad the actual geometrical dimensions of the core are used. ‘The signi- 
ficant dimensions and coontinava system used are shown in Figure 1. The 
numérator of haraeereeston for the reactivity contains terms which repre- 
sent changes in particular reactor parameters in the perturbed system. 
The perturbation is assumed to occur only in one disk of the core at a 
time and is uniform amin erat disk. The total reactivity of thet déte 
is found by summing the reactivities of oeas disks. Since the perturba- 
tion is zero. everywhere in the core outside the disk under consideration, 
the limits of integration in the numerator are from zero to R in the 
| yadial direction and from h to h + Sh in the vertical direction. The 
denaaidacor,: noweveds contains no perturbation terns and is integrated 
over the entire volume, i.e., from zero to R and from zero to H. 


ae | Le 4, SP VECE, B* AV 
y= Sy scvers, p*av 
1,7 J) 38, verc o* dv 
Lye Sy 88, 0%AV 
Lge J adver av. 
T= S) vercdstav 


Therefore, : 5. bs Ti+ I.+ 4 fi Tu -Is = : (19) 
b 
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Fig. 1 


Fuel Core Dimensions 
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Evaluation of these I terms is presented in Appendix C. Substituting the 
resulting expressions for the I terms in (19), and cancelling like terms 
results in 


§- SPVECEL + SCVEPEF + SE/VEPCQ ~ $E0Q- sD [o*aq+(F)y] 
TVEPC Ze 


This expression can be further reduced by factoring out the $ term and 
dividing each term in che: aunerator by the denominator. This step is 
quseieied only if it is assumed that the values of P, C, 2, %_, and D° 
appearing in the numerator, which refer to disk properties, do not vary 
significantly from their corresponding values in the denominator which 
refer to the entire fuel core. Making this assumption makes the final 


reactivity equation much more tractable. Then, 


3° $0 % 320 MAP 
5" 21 %+ c +} . Deets, DEPOSs [s*+(F) ]s} 


F 





vou : Ze 
Since K,* HEPF where f* and %* Yar the expression for § 


can be written as 


Qsye . Ce . Be Fe [ zr 
Ss i T fp 5 C i 2s » Ree "Ze mon “ or £490] 
Toyz . 
Letting, e2 X*4 Cry 7 


and since Lees Za? the. final expression for the reactivity for this 


case is 





+ 


= 2(# So, SB _ 1 Ado _ Le 2 ) 
Le ae, ee eee, Sa Ree 
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2. Extrapolated Dimensions 


The derivation of the equation for the aeticiey for the case 
using extrapolated dimensions is similar to that of the bare core geo- 
metry. The difference results from the fact that for this analysis the 
flux is no longer considered to be zero at the physical boundary of the 
fuel core but instead goes to zero at a hypothetical distance, d, beyond 
the actual boundary. The resulting flux distribution more accurately 
represents the actual situation within the reactor core since experimental 
evidence indicates that the flux is not zero at the core boundary but has 
a finite value. It can be shown by transport theory that for a planar 


free surface, d has the value 


A 2 0.71 Ny 


The transport mean free path, _ Nar -is found from the equation 


As 


Me? 
sue tr }- aH 
But since 4#{ is small for the polyethylene moderator, it is assumed 
hat er \s Ls 


Substituting this into the equation for d gives 


0.71 
as 


For a curved free surface, such as the outer edge of the cylindrical 


d= 


fuel core, the above equation for d can still be used and introduces 
negligible error provided that d is small compared to the radius R. 
This condition will be shown to be true for this reactor. Therefore, 
because of the extrapolation distance, the extrapolated radius and extra- 
polated height of the fuel core become 

R° =R+ d 

H’ =H + 2d 


and the parameters % and ® in the flux expression become 
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From this point on, the development of the reactivity equation follows that 
outlined in the preceding section with two exceptions. First, several 
‘of the zero order Bessel functions in the reactivity equation for the 
bare core geomet ry case were zero because of the boundary condition that 
the flux equal zero at R. When using the extrapolated radius R' these, 
Bessel functions are not equal to zero and are retained in the equation. 
‘Secondly, using the extrapolated height H', rather than H in many of the 
facegvale: tavolving trigonometric functions yields a slightly different: 
result. Because much of the development of the reactivity equation is 
repetition of the bare core geometry case it will not be repeated here, 


but it canbe shown that for this case that 


Te, 30, ee | a _ th 
ATR EN Ae ose i. ZA (= . GG . Zp Keg Zea SS Kees ) 
where . 
. Uh eo aitCneSny Let, 2th 
Te — sh = => sin = 
TTH «7th 
A= arte sin = 
2 ae ty we 
be me Gi) F 
= Tan 4b gin STH 1 sin AUChtSh) 
We : r + Zz Sin 4 z Sin i 
K* J (aR) - J, CKR) TJ, (R) 
aaa E* J,*(«R) + J,* CR) 
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3. Reflected Core 

Since the fuel core of the AGN-201 reactor is surrounded by a 
graphite reflector, the expression derived for the reactivity in this 
section probably is the most realistic and accurate of the three cases — 
discussed, even if it is only an approximate selacieae. As it. turns out, 
the problem of a cylindrical core with radial and end reflectors, as is 
the case for this reactor, is more than can be handled in an exact. analyti- 
cal manner. The problem arises from attempting to find a form in which. 
to represent the flux distribution so.as to satisfy the boundary emda 
tions at the “corners” of the cylinder, i.e., at r = R and z = H/2. The . 
problem can be solved by numerical methods using réelexation techniques 
but is wether complicated and will not be used here. 

An approximate method of solution is to treat the uaceor as having 
a radial reflector and bare ends but having a height of H + 2% where 3, 
is the end reflector savings given by : 

$,2 -@w tan (SE tonn <4) 

Therefore in the flux distribution expression, 


an 
A= H+ 232 


Since the reactor core is considered to have only a radial reflector 


os 2:408 
2 R+ Sr She . 
— Spe + tan”! (FeShe tanh = 


From Appendix A it can be noted that the thickness of the reflector is 


uniform about the reactor core so that 


1 
From this, it follows that 
$,% 3,2 3 


Following this preliminary discussion, the development of the reactivity 
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expression follows closely that of the extrapolated dimensions case, 
only consideration is directed toward the effect of the reflector savings 
concept rather than the extrapolation distance. 


The resulting expression 


for the reactivity for this case is 


S= (3+.  - | $24 _ .L*a 80° 
2G\P  ¢ Zp Ke Za Kot 9 





_ Th 2 4 es, 2MCH+$h) 1 - 2ity 
nee PS Gog En Wes za 
_ 1TH os. OTH 
G* saat a8!" Tees 
— xa® wT 2 § 
ae E Cre F 
and 
Pee A So «jy, 2 2 1 amCne$h) 
S H+as oe 20 Heads zn H+2$ 
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IV. COMPUTER ANALYSIS 


A. DERIVATION OF NUMERICAL COEFFICIENTS 

In the preceding sections, equations for the change in reactivity 
of the perturbed reactor system for three different models were devsioned: 
These equations were: 3 


Bare core geometry, 
- 2/2, Soe ie ye 5e) 
fone S. KeoC fa Kec D _ @ 


Extrapolated dimensions, 





a Te, 00 ee LL She | th 90, 
5 an (+e te Ke® “fq ~ Kao =) (24) 


ge E(B 3, ew) 
S* z6(% CS, KO 2a. Kee x) ee) 


As can be seen, each of the equations contains the. terms P, C, Zes Kw 5 
z faa and D. These parameters citer to the cylindrical fuel core 

as a whole, aré constants, and are evaluated at their initial values, 
prior to the reactivity drop.. From Appendix D, where sepranalone for 
_these parameters are developed and nmnerieatiy evaluated using data 


from the AGN-201, the following results are obtained. 


P = 0.9473 
C = 0.7414 
Z_ = 0.06507 om” 
Kw = 1.510 
= 0.1010 em” 
to BAIT ea 
D = 0.2199 em 
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Before these reactivity equations describing the three models can 
be used, a discussion of the term $C must be presented. Recalling that 
Ze “178 
c= = ta = ) 
@ m (=, : 
where B is a constant and 
$ = files 
eo” 6 GA) 
the problem is to determine the manner in which t, ‘varies as a func- 


tion of the moderator density. The Fermi Age, - » is proportional to 


the inverse of the product of the logarithmic energy loss per collision, 





& » and the macroscopic scattering cross-section, Zs ‘ 
| 
@ ae anemia 
® 
5 2s 
To make this relation an equality, let the constant of proportionality 


be €, , so that 
€, 


| er 
From Appendix D, C = 8.459 eae 
| # = 0.9090 
and Zs = 1.518 em” 
Therefore, €, = 17.75 


Since _. 2, is a function of the moderator density, the relation between 
e 
cS and hence C, and J" is. 
. 17.75 
: Fzs 


1. Bare Cy lindriéal Geometry 


The term e in the bare core geometry equation (20) has been de- 


fined as 

e=a'+(H) & 
From Appendix D, o* = 0.03587 cm 
and H = 12.7 cm 
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Therefore © = 0.03587 + 0.1695 “@& (23) 
Substitution of the values for the core parameters and (23) into equa- 
tion(20) yields the numerical equation for the reactivity for the bare 


core geometry model, which is 
$= D[ci.056) $P + (1.899) 8C + (18. 37) 9.2 - (8,849) 9S. 


~ (0.3172 + 0.1497.) SD] £24) 


2.  Extrapolated. Dimensions 


In the equation for reactivity in the extrapolated dimensions 
model, (21), b represents the quantity 
oa xo* ott \* w , ake 
be Ge (25) 


From previous analysis it has been assumed that d, the extrapolation 


distance, can be calculated from the equation 


On 


| a 
From Appendix D, 2s # 1.516 ene 
so that = | d = 0.47 om 


The extrapolated dimensions are therefore 
R' =R +d = 13.17 om | 
H' = H+ 2d = 25.07 om 
In equation (25), for the extrapolated dimensions case 


2.408 = 
er = 0.1826 om 


J, GR) = J, CHR) J, CR) 
Jo* C&R) + J,* C&R) 


1 


X 


and E 


From a table of Bessel functions, and since & R = 2,319, 


J,4R) = 0.04531 


J, C&R) = 0.5364 
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J,@R) = 0.4173 


So that X = 0.5175 


and | E = 0.2898 | : 


Substituting these values into (25) yields 
, Ww 


b = 0.05954 + 0.01631 = (26) - 


Inserting the values for the core parameters and the expression for b, 
(26); into equation (21), gives the numerical equation for the re- 


activity in the extrapolated dimensions model as 





$= Fp [(rose) se + c1.349)8C + (16.37) 82, - CB.944) $Lq 


~ (0.5265 +0,1442 ) 5p] (27) 


3. Reflected Core 
To evaluate the numerical coefficients in the reactivity 
equation for the reflected core model, (22), the first step is to ve 


evaluate the reflector savings term, % . 


iss ste ite eee gS Se oy 


Q= + ten” C Ober tanh +) (28) 





From Appendix A, T = 20 cm 
The values for B and D have already been given as 

B = 0.2298 cm”! 
D = 0.2199 cm 


Since the reflector material is high grade graphite, 


L. = 54.2 em 
r ‘ i« 
D. = 0.94 cm 
r 
Therefore, evaluation of (26) yields : 
% = 3.49 cm 
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so that R + 3 = 16.19 cm 
H+2% = 31.11 cm 
For the reflected case, the parameter Q& in the expression for the flux 


distribution has been given as 


1 





2.408 i 
= Gos = 0.1485 cm 
so that a %R = 1.887 
Since » ~ x = a (R)- J, CAR) I, (aR) 
cand _ B® T2@e)+ J 2(ar) 


“Substitution of the values for the Bessel functions of 
Jo@R) = 0.2898 
J,@R) = 0.5814 
| J,@R)= 0.3269 
yields _ X = 0.4868 
E = 0.4218 
. Substitution in the expression 
| | ket 7 ES 
* as “e+ (qs) F 
gives an expression for a of 


Ss seen 5 
a = 0.02597 + 0.01020 a : (29) 


woe 
e 


Putting this expression, (29), and the values of the core parameters 


into (22) gives 


S= Fe [(.060) $0 + (1.349)8C + 15,37) 8p - (8.944) BZ 


- (0,229%+ 0.09620 $) 30] Go) 


as the numerical equation for the reactivity in the reflected core model. 
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B. . FLUX WEIGHTING OF CHANGE IN MODERATOR 
It will be noted that four-figure accuracy is maintained in the 
ava laation of all the above parameters, This is not strictly justified es x 
on the basis of the accuracy of the data available on the AGN-201 re- 
. actor, but was adda weheaaey by the very. small magnitude of the reactiv- 
ity change, 0.00073, which was to be accounted for. The validity of the 
above numerical results is strengthened somewhat by using a weighted dis- 
tribution of the change in density within the fuel core. It seems al- 
most certain that the change in hydrogen density is not uniform within 
the fuel core, but is instead higher at the center, where the neutron 
flux is highest and decreases towards the edges of the core. While the 
density change was assumed uniform within a single disk within the core, 
the distribution of the density change in the wereieel direction was 
"weighted" by a cosine function corresponding to the shape of the flux 
distribution in the vertical direction. 
Define a weighting function = cos Ge 
Using the distance of the centerline of each disk (24) from the zero 
ordinate of the z-axis of the fuel core, the fraction of the total 


density change within a single disk Ag » can be found from the expres- 


sion i 

005 “iT ee e 
Ba Te a ean 
Zs C05 4 A" 


na 


The results of applying this equation to the data listed in column two 


of Table 1 is shown in column four of that same table. 














TABLE 1 





Disk 24 cos Lal ; A 

1 | 4,66 0.02973 0.00584 
2 4437 0.1253 = 002464 
3 | 4.00 0.2455 | | 0.04828 
4 3.44 0.4198 0.08256 
5 2.68 0.6323 0.1243 
6 1.94 0.8012 0.1576 
7 | 0.83 0.9626 _ 0.1893 
a ~1.00 | 0.9458 0.1860 

9 “2.50 0.6773 0.1332 
10 4.00 0.2455 0.04828 


C. COMPUTER PROGRAMS 

Al11 of the discussion in this section to this point has served to 
establish the basis for the final application of perturbation theory. to 
the analysis of the behavior of the AGN-201 reactor. Appendix E con- 
tains a digital computer program, the results of which demonstrate the 


reactivity change of the reactor as a function of moderator density change 





caused by operation of the reactor. In ihe: peontaan the change in moder~ 
ator mass starts at 10 pm and increases to 50 gm, in increments of 10 gm. 
This overall change in mass is weighted according to the results listed 
in Table 1, to find a dewatty change per disk. For each change in 
moderator mass, the program uses the denceistey-watakous (24), (27), 


and (30) to calculate the reactivity drop per disk and for the entire core 


4i 








for each of the three models describing the reactor core. The numerical 
results of this program are presented formally in the next section. 
However, at this point, it is sufficient to state that the density ca 
change predicted by the perturbation theory equations corresponding to 
the actual value of the reactivity drop equal to 0.073% are as follows: 
Bare core geometry: AM, = 17.9 gm 
Extrapolated dimensions: AMp)= 15.1 gms 
Reflected core; AM, = 13.8 gm 
These Sanger tes the total core mass change are used in the computer pro- 
gram presented in Appendix F, where the effect of adding a new fuel disk 
of known properties is analyzed. Using the data for the new fuel disk 
listed in Appendix A and the results of the first computer program to 
determine the present properties of the old disks currently in picags 
actor, the resulting excess reactivity of the reactor is diated: dink 
the new disk displaces a disk of similar size by again using the numeri’ 
cal perturbation equations (24), On, and (30). The results of this 
vist te, t 


program are presented and discussed in the following section. _ | i 
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V. RESULTS AND CONCLUSIONS 


The numerical eaiiis of Computer Program 1 presented in Appendix E 
are listed in Appendix G and: plotted graphically in Figure 2. te ap- 
peace from ‘Figure 2 fier. a linear relationship exists between the reactiv- 
ity « of cle! ‘reactor core and the change in moderator sdaapity within the 
cores, Taking the slopes of these plots for each of the three models shown, 
‘the: following espreestons can. be derived:.. 

Bare. core geometry: AS = -0.00004095 My 
Extrapolated dimensions: AS acco: 00004888 Me 

Sie . Reflected core: ‘Ade ~=0. 00005302 M 
Kivoy: from Figure 2.it is possible’ to obtain the change in polyethylene 
ass 1 for. “Sech model which corresponds to. the cpeerves change in reactivs 
eee of the reactor. of 0. 00073. These: changes are. ; 

ae core Se: = 17.9 gms ja 
mittapotaced: dimensions: M, = 15.1 gms 
| Reflected core: Mo . 13.8 gms 
In previous discussion a value of the weight of p hydrogen bled off from’ 
‘the reactor during its lifetime of operation: was” ‘estimated as 23.14 gm: 
While this value. is slightly ee ae that predicted by ‘the perturba- 
tion, theory equations, no true ‘error analysis can be conducted. The 
value of 23. ” is actually only | a crude approximation in ‘itself, but the 
fact that this -experinental value and the theoretical palien are reason- 
ably close to one another lends. validity to che pewutes of this analysis. 
. It can also be stated that the savorat® erux level of 20 watts used to cal- 
culate the value of 23. 14 gm must be alone. to the actual value of the 


reactor flux. Since this is the case, bine arguments discounting fuel 
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burnup and reactor poisons as sources of reactivity drop are reinforced. 
Interpolation of the data in. Tables 6 and 7 of Appendix G to corres- 

. pond to the above changes in moderator mass allows an analysis to be 
made of iheHenctawies drop in the reactor core as a function of vertical 
position. These results. are procrce in Figure 3. A6 ‘can be seek, the 
curves for the three models are different. “The ‘extrapolated dimensions 
seer curve, in general, siieard a | higher overall value for the reactivity 
decrease than the bare core model curve. The difference: ‘noted between® 
the reflected core model curve and the extrapolated dimensions model 
curve is that the reflected core curve is slightly “flattened”, The 
reactivity feecaias is higher at. the top and bottom of the core and lower 
in the middle for the reflected core. The variations among the three ~ 
curves reflect the various shapes of the assumed flux within the core “~~ 
for each model. 

‘Using the density changes for the three models of 17.9, 15.1, and 
13.8 gm to determine the present core density of the reactor, the re- 
sults of Computer Program 2 in Appendix F show that the effect of inte#s’ 
changing the new disk with an old disk of the same size, (either disk 
4,5, or 6), results in a emai! increase witch exceeds the ligensed 
limiting value of 0.40%. ‘This can be: confirmed by referring to Table 2 
or Appendix H. The 1acee reactivity increase can be attributed to the 
large increase in moderator aenaley Hekeaen the old diske ‘and the new 
disk. This increase in density is approximately 0.2 gm/om” higher for 


the new disk. By modifying Computer Program 2 to print out the values 





of $C and %D, it was found that these two positive terms in the 


reactivity equation, (§D is negative, but is multiplied by a negative 





coefficient) which are functions of moderator density, increase substantially 








Reflected Core 





Model 

-0.0015 
a 
n 
a 
S 
a 
oO 
= 
BH 
ae Extrapolated 
a Dimensions 
A Model 
td 
z -0.0010 
i] 
vu 
a 
= 
Nw 2S 
6 Bare Core Model 
% 
a 
u 
Qo. : 
A» -0.0005 ° 
> 
Ww 
A 
> 
ord 
Ps 
Q 
g 
pa 

0 





4 2 0 a) 4 
Distance from Core Centerline 


Fig. 3 


Spatial Variation of Reactivity Within Core 


46 








faced Sie cia anes gate: 
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REACTIVITY 


Extrapolated 
Dimensions 


0.004392 | 0.004816 


0.005938 


1. 0.006481 0.007351 





TABLE 2 
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0.006643. | 


| 0.007707 


Reflected 
Core 


0.005767 














when the new disk displaces an old disk. These increases more than com- 
pensate for the negative contribution to reactivity caused by the de- 
crease in fuel density of the new disk as compared te the old disks. 

For purposes of illustration, consider the values obtained from Computer 
Program 2 listed in Table 3. These values are those obtained if the new 


disk were to displace Disk 4 using the reflected core model. 


TABLE 3 
$P = 0.006269 
@C = 0.04459 
SO, = 4.931 b 


$B, =-0,.00003234 cm”! 
Soi" = 5.892 b 
300s 0.02357 b 


326 = 0.004563 em! 


$%_ = 0.004523 em” 

20 =-0.02916 cm 
From this table it can be seen that the changes in parameters caused 
directly by moderator density. changes, when multiplied by their corres- 
ponding coefficients in che ceactiotes equation, are of about the same 
order of magnitude as those caused indirectly by the change in cross-~- 
section caused by the change in neutron temperature discussed earlier. 
Thus, it is concluded that a simple interchange of the new fuel disk with 
any one of three similar size disks which are presently in the reactor 


core will produce too large a value for the excess reactivity. 


‘The conclusions drawn from this perturbation theory analysis of 


reactivity may be itemized as follows: 











1) The decrease in reactivity due to loss of moderator within 
‘the core is a linear function of the change in core density. 

2) Due to operation of the AGN-201 reactor over a period of 
several years, the moderator mass has decreased from between 14 ‘to 18 -gms, 
a 3) The interchange of the new fuel disk with an old disk .of 

similar size currently in the reactor and in its present position within 


the core will produce an excess reactivity exceeding a value of 0.407%. 
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se 


VI. RECOMMENDATIONS 


- SRST 


“There are several areas into which further study of the problem of 
the reactivity analysis of the AGN-201 reactor could be directed. Based 
on the results obtained in this thesis, a direct displacement of an old 
disk of the same size by the new disk would exceed the limit of 0.40%, ex- 
cess reactivity. The perturbation equations developed here could be used, 
however, to predict the change in reactivity resulting from any new combina- 
tion of disks in the core, such as displacing an old disk with the new 
disk but removing the smallest top disk, adding the new disk at the top 
of the core, etc. Only slight modification of the data input into the” 
computer programs would be necessary to permit this. In other words, the 
equations sila conndter programs developed here could be deed to find a 
disk arrangement. such that the excess reactivity approached 0.40%. 
Secondly, a chemical analysis eéuld be conducted both to determine . 
the exact products formed after polyethylene is irradiated in the core’ 
and the nature of products present: in the evolved gases. In this thesis 
the assumption was made that only hydrogen was evolved gadvthar no new 
type of molecule was formed in the fuel core. Perhaps a chemical analy- 
sis of irradiated polyethylene could either substantiate this or illumin- 
ate more clearly the exact nature of this dcataanen 
Third, the method used to relate neutron temperature to moderator 
density, equation (8), was only an empirical relation which was modified 
to suit the data in this thesis. A more exact method of handling ‘this poof 
problem is to treat the problem by an energy-dependent model, with the 


scattering kernel computed by the code GAKER, which is available through 


such agencies as the Argonne National Laboratory. The use of such a code, 








if obtained and programmed to run on the IBM-360 computer would undoubted- 
al ly increase the accuracy of the results of this type of analysis. 

: Finally, other analytical improvements could be made in this analy- 
sis. The results of a chemical analysis, such as previously described, 
might indicate that a two-group analysis rather than one-group as used 
here, might provide better results. .Also, a more exact analytical method 


of handling the reflected core model would be desirable. 
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APPENDIX A 
DESCRIPTION OF THE REACTOR 


A. | GENERAL FEATURES: 

_ .The AGN-201 reactor is a low power nuclear reactor designed pri- 
marily for research and educational purposes. .It. has a fuel core con- 

. sisting of ten uranium impregnated polyethylene circular disks 10 inches 
in diameter sna gargtiie in thickness from 3/16 inches to 1-1/2 inches, 
stacked on top of one another to form a cylindrical core. The uranium 

-is in the form of uranium dioxide enriched to 19,937 with U-235 and is 
uniformly dispersed in the polyethylene. Inserted in the core through 
its base are four control rods. Three are the same size; the earth is 
smaller. Two of the larger rods serve as safety rods, the third as a ~~ 
coarse control rod, and the smaller rod functions as a fine control rod. 
Since these rods are of the same material as the core, their contributions 
to reactivity considerations are handled by incorporating their properties 
into the disk of which they form a part when fully inserted into the core. 
This is both justified and necessary when calculating the reactivity in- 
crease caused by adding a new fuel disk since insertion of the rods in- 
creases reactivity and the maximum excess reactivity value is not to be 
exceeded, 

Surrounding the fuel core on all sides is a reflector of high density 
(1.7 gustan) graphite, 20 cm in thickness. Shielding is provided by a 
10 cm thick lead shield, a thick steel tank having a radius of 43.4 cm, 


a borated-water-filled shield tank 6-1/2 feet in diameter, and finally, 


a housing of large concrete blocks which encloses the entire reactor to 








provide protection against radiation while operating at higher power 
levels. The general features of the AGN-201 reactor are shown in Figure 
4. The core tank of the reactor is fluid tight and leak tested.to a seeas 
sure of 5 psig. To handle pressure build-ups divine Seetation wad after 
shut-down, the tank is connected with copper tubing to a gas hand tine 
system which provides a means for measuring core pressure and a means for 
bleeding, storing, and monitoring the gases generated by the fuel disks, 
To insure against ignition of the hydrogen gas with any oxygen in the air 


which might have been initially in the fuel core, a system feeding nitro-. 


gen, an inert gas, into the core tank is provided. 


B. THE FUEL CORE: 

The dimensions of the fuel core bia of the ten disks composing it 
are shown schematically in figie S. For conven sues of reference, the 
disks have been arbitrarily assigned numbers. Pertinent disk and core 
data is listed in Table 4. It should be noted that while the various 
physical ppapeciies are constant within a single disk, variations exist 


from disk to. disk. 


Cc. NEW FUEL DISK 
Data pertaining to the new fuel disk to be inserted in the fuel core 


of the. reactor is presented in Table 5. 
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TABLE 4 


ha reacting 


Fuel Core and Disk Data 














ee, ee hag 7235 | Pigasgl 
1 ss -292.5. 224.6 1.196 0.06958 0.2796 : 
2 — 498.2 482.7 1.032 0.05625 0.2449 
3 552.4 482.7 . 1.144 0.06648 0.2670 
4 1038.1 965.4 1.075 0.06259 0.2515 
5 1049.1 965.4 1.087 0.06313 0.2536 
6 1047.2 965.4 = 1.085 0.06303 042532. 
7 1721.2 1930.8 0.8914 0.05179 0.2080 
8 1722.9 1930.8. 0.8923 0.05190 0.2084 
9 1761.8 1930.8 0.9125 —«-0.05298 0.2129 . 
100«1784.8 | 1930.8 0.9244 0.05371... 0.2157... 
Total core weight = 14821 em 2 | : 
: | 


Total core volume = 12230 cm 
Overall core densities: 


PP = 0.9383 gm-cmn™> 


F 235 3 


= 0.05451 gm=cm™ 


POO 26 -i0d: gneen 


S26 19736 easen © 


TABLE 5 
New Fuel Disk Data 


‘Disk volume = 965.4 gen 
Mass of polyethylene = 921.4 gm ‘s 
' Density of polyethylene = 1.268 ‘gm-cm 
Density of U-235 = 0.06003 gm=cm., 
Density of U-238 = 0.2287 gm-cm 


‘eon 
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APPENDIX B 
DISCUSSION OF THE ADJOINT OPERATOR 
The fuel core of the AGN-201 reactor is a finite cylinder, and the 
diffusion equation governing its. behavior is a. partial differential 
equation. While the following discussion applies equally well to such 
a case, for simplicity of presentation, consider the following differ- 


ential equation which is the one -group diffusion equation for a critical 


slab reactor shown below. 


A ney, oS i 
oR BOO GR + FCA) ® O 





This can be written as 

MG=0 GD 
where M is the operator 
eine ie M = ax DX) ay * FO) 
The equation ae oe 
is related to equation (31) and is referred to as the adjoint equation. | 


M+ is known ag the adjoint of M and the function # is referred to as 





the importance function. In addition to satisfying equation (32)» y 
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must satisfy the same boundary conditions as @ , which for this case 
is to vanish at x = * %z » the slab half-thickness. The adjoint opera~ 
tor M+ is defined as that operator which satisfies the equality, 


G/z, Oy 
f vMw ax> Jy wM*y ax (33) 


~My - 


where v and w are any two functions which vanish at X = tM . Sub- 
stituting the expression for M into the left-hand side of (33) and 


integrating by parts yields 


ap. a a . : j Oy, 
ie vo ax Dax *E]wadx = CvDw’-v x) a 


Ore 
eles w La D a *Flyvax (34) 


te 
But since ; v¢t Ys ) zo 
and ‘ W (? 2) = Oo 
; ‘ “Va. = 
then, ' CvDw'- Vv Dw) hay = 0 


and from (34) . 
2. Oye, Fe 
J vMwdx = ee whvadax (35) 


Comparison of (35) with (33) shows that for one-group theory 


M+ = M 
When an operator and its adjoint are identical, the operator is termed 
self-adjoint. Further, because equations (31) and (32) are homogeneous, 
and F and D satisfy the same boundary conditions, @ and ¥ 


are everywhere proportional. [5] 


58 








APPENDIX C 





EVALUATION OF INTEGRAL TERMS IN REACTIVITY EQUATION 


Evaluation of the I terms in the reactivity equation (19) is as 
follows: 
L=S) spvecz, gt av 
For a cylinder, av = 21 cae aL 


and the flux distribution, D >» is 


p= Do a (ar) aos BL 


where 
2.408 
| ae ae 
1 
and * is 
Therefore : 


ne Sh 
LJ te sPVECS, 9: 38°) 00s 2 (ezXanrde AZ) 
Since UV, €& , C, Bes arid ®. are constants within a disk and SP, 
the perturbed parameter, is a function only of a density change which 


ig assumed uniform within a single disk, 1, can be written as 


en 
I,* 2 sevecd, BES cIg (andar fr oos* (62) AZ 


Evaluating, 
. 2 6 Le 
Se rtte@rac® £I%@R) 
since J,* CaR) s a ( aHg6 8 ) =O 
ae HQ 
J, eos (62) AZ = ar. 
where - Tt ; an (h+Sh 0s 
: : | qs + a sin aH -t sin 
merefore, g 1, * % WPVECA PRT R) HO 


Through a similar analysis it is found that . 
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1,2 Z 2CVEPE, BERT, CAR) HO 


I, = & 32, VEPCQ R*T,*(@R) HQ 
aa Ty %23q G2 R*T,* CAR) HQ 
The solution for the I. term is somewhat different since the diffusion 
coefficient is weighted by the gradient of the flux rather than by 
the square of the flux. \ The gradient of the flux can be written as 
WU, Bra tse +u,¥ 
But, ty F 4 =O » Since it is assumed GF Z(O) 
Therefore, . Cve)* - ()° . (3¢)* 
Taking the derivatives of p with respect es r and z, and squaring 
yields : 
(VG) = GP [8 Tr) cos* (62) + B*J,(#r) sin* (4Z)] 
Therefore, I, can be written as 


R 


is* am300. {J 


+3 
Bt {ete cryar +f 'sin* (62) 02} 


ee nth ; 
rJPcarddeS, cos*( 42) Az 


a R a emanation Ry" a 
% i rJ, (andar * = J *(#R) 
heth ‘ HO 
a. cos *(@Z)aZl* am 


R ee 
Jo cd, Caedae = 2% T*(4R) 


neh o HY 
Je sintcozvar= & 
Th ath . t . amCn+sh) 
where You z sin ~ HO” 2 Sin 4 
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Therefore, 


(AR)> 2 
Lee zmsoge § S° 5,*caR) SE + GH) Z T2 CR) 


oe R phat 2 rn 
Ips ¥BDGZR*T CR) H LV Q+ CH) YJ 
The evaluation of the integral in the denomenator of the reactivity 


equation is . 
ie J, vere 3,G°AV= ZTVEPCS, BF [Err°@ Rar te cos COZ)dZ 


i cd, (ror = © 5% (ae) 


ud 
Sf" cos* (#2) az = 


Therefore, . 


I= 2 vePC2, 62 R*T,*(#R) H 


61 





APPENDIX D 
NUMERICAL EVALUATION OF REACTOR CORE PARAMETERS 


A. RESONANCE ESCAPE PROBABILITY 


A general expression given by Glasstone [4] for P, the resonance 


escape probability, for a thermal-neutron reactor is 


Za SE. rs 
Pca) exp L¥ ft 4 e (36) 


For a homogeneous reactor system, in which Za is assumed to refer to 
the fuel only, and Z. is taken to be independent of neutron energy in 


the resonance region, (36) becomes 


Eo ae a 
PCE) * ™ kp [- ¥, Fy, i (ark zs o “) = a! (37) 


The integral in (37) is called the effective resonance integral and 


I: ns ( BE, % q.” ca 


Rewriting equation (37) 


can be represented by 


eee? 
P= exp [- os =] (38) 


For the AGN-201 reactor, 
P uv e.FP vg 
+ = 
2,*2, =NG +Noa, 


: P u 
Dividing by N, 


NP 
Fe Teale 
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v = 
However, since J, is small and can be neglected with respect to the 


first term, 


< NP fp sf 
Ne Ne EF AH (39) 
ef 
P' Na =e 
. 6 ss 
Because es “BP Ts 
and NY = ao 
equation (39) can be expressed as 
Zs 5° 7h Pye is 
Spt ar (A , (40) 


From Appendix A, 
e Ss 

SP 05679383 guicem”> 
: we ~: 

and 3° «0.2736 gm-cn™? 


&. e.. , 
The evaluation. of (= ) is as follows. The exact molecular weight of 


_ polyethylene, aP is unknown but can be written in terms of the (CH,) 
groups forming the molecule as 
AP = na# © 14.0 n (1) 
Similarly, 3 can be written as . 
; 7 
T*? NG" + ZN, (42) 


Dividing (42) by (41) eliminates n and hence the exact molecular weight 
a 

need not be known to find ( a yP. This same method is employed to 
find the ratio of absorption cross-section to molecular weight in a 
subsequent paragraph. Correcting the cross-sections in (42) to an 
effective neutron temperature of 340° K through the relation 

gt. ini ye 

U3 * TIZ3 V tee Ts 
yields 

@° = 3.35 b 


t"* 16.8 b 








“= ne 24° 
(By. Bee eer 


Therefore, 


Also, i 
A = 270 


Substitution of these values into (40) yields 


“Ne 7 2472 (43) - 
An emperical relationship for the effective resonance integral in which 


U-238 is the resonance absorber is given by Lamarsh [5] as 


I = 2.73 ( (44) 


23 ) 0.486 
N° 
Substituting (43) into (44) 

I = 121.7 (45) 
The average logarithmic energy decrement per collision, é » fora 


mixture of compounds such as in the fuel core of the AGN-201 is given 


as 


- Y ¢ « 
§ = aS a 48) - 


v ¥ e é 
Zig * or = 0.0084 cm” 
where . a = 13.7 b. 
VW 
§° = 0.036 


e e, ¢ SaP Z 
By = S°NaCH) = 1.508 em? 


The expression £08 P is — ec ws 
; ef, (8%) + 2066) 
g, ¢ + 2," 
c 
Using 5 = 0.158 





4 
§" = 1.00 
the value for § P is 


gf = 0.914 
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Substituting these values into (44) yields 


EF = 0.9090 | (47) 
Inserting (43), (45), and (47) into (38), the value of the resonance 


escape probability term is found to be- 


P = 0.9473 


Be FISSION CROSS-SECTION 


The equation used to calculate the fission cross-section, Z es 
236 


; ase Na_ x 

is | =, = pas SF 
a3g 23 

From Appendix A, Sr" 0.05451 gmecm 


The microscopic fission cross-section, corrected to a neutron tempera- 


ture of 340° K is 


- Gs 466 b. 
Since N= 0,602 
235 


and AT" & 235 
the fission cross-section for the overall core is 


3. = 0.06507 ca 


GC. INFINITE MULTIPLICATION FACTOR 
The expression for the infinite multiplacation factor, K,, , is 


known as the four factor formula and is given as 


Kg ere 
The value of % can be found from the equation 


$, : 
m2 pe? St 


a 
Since in this expression, the “fuel” refers to the mixture of U-235 and 


U-238, Ba is given by 


e. ass ass 
Zia” 2 = de 











235 : 
Therefore, Aen se gm 


Since py 235 2.47 


238 = 258 . 7 
379 Jone G = 0.07764 cm 


1 


_ 235 
where Ta = . 556 b 


23% 
Nea =- 
and eee ee 


233 sa']. 
“yar Me = 0.001236 cm 


4 = 2.038 
It is assumed that the fast fission factor, E , is equal to unity, 
based on the assumption that fast fission of U-235 does not contribute 
to the multiplication. From a previous paragraph 

P = 0.9473 


Finally, the equation for the thermal utilization term, f, is 


f , S32 , _S0 : 
Zo soe at 
2 t8S | (288 ES 
3. a+ Za, = 0.07888 cm 7 
ay : 
sf = §*Ng (SP) = 0.02209 em™ 
: Ge yf 
where (=) = 0.0391 b. 


The resulting value for f is then 


£ = 0.7810. 
Inserting these values into the four factor formula for K,, , the 
infinite multiplication is found to be 


Keo = 1.508 
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D. ABSORPTION CROSS-SECTION 
The total absorption cross-section for the fuel core is found by 
summing the cross-section of all the substances present in the core. 


Therefore, 


6 


oes ee eat Cog 
a . Zs + o ag me ) Za x Da 


Sa = 0.1010 cm”? 


E. DIFFUSION COEFFICIENT 
From an earlier discussion concerning extrapolation distance, the 
assumption was made that without introducing significant error to the 


calculations, D could be defined as 





Gg 


og ieee tee Me ei be 7 Best 
where 2g * 2. * Zs, = “ATG, * §No (+) 
“since —, Gy = 13.7 b, a 


3,.=- 1.516 em” 


and the diffusion coefficient equals 


D = 0.2199 cm 


F, DIFFUSION LENGTH 
From diffusion theory, the quantity referred to as the diffusion 
length is defined as_ 


vos 





B 
. Za 
Since Za and D have both been previously calculated 


w= 2,177 em” 
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G, FERMI AGE TERM 
The Fermi Age term must be evaluated last because it is solved in 
terms of all the other parameters in the criticality equation. Recall 
that the quantity C has been defined as 
Ze ~I é : 
c= -8 ton! Cz) (48) 


For a cylindrical reactor core, B is given by the relation 


pe = t+ 6° 
ZMO08 -1 
where aA —_ 0.1894 cm 
and 6 = i = 0.1302 cm”! 
Therefore, B = 0.2298 
By definition, 
Z. = aly 


The problem involved when dealing with an hydrogenous moderator, 
such as polyethylene, is that Le cannot be calculated exactly from 
the Fermi Age equation and approximate methods must be resorted to. If 


one were to write the critical reactor equation 


KeG 5 | 
| + L*6* 
to include the built-in excess reactivity of 0.000354 of the core, the 


resulting "supercritical" equation would be 


se = 1.0004 


Since the only unknown in this equation is C, the solution is 
C = 0.1704 
From (48), 
wt 
as tan (+) = 0.1704. 


Since B is known, this is an implicit function in terms of e. - An 


iterative solution yields 


68 








at 
Zo = 0.1985 em 


r> 


hi 


Therefore from (47), 8.459 et 


Q 


H. SUMMARY 


The results of the preceding paragraphs may be tabulated as 


P = 0.9473 
2. = 0.06507 a 
Ko = 1.508. 
a? 0.1010 em”? 
D = 0.2199 cm 
wv = 2,177 aa 
C = 0.1704 


se 


. 8.459 er 








APPENDIX E 
COMPUTER PROGRAM 1 


COMPUTER NOTATION 

AS = owe neutron temperature of disk 

A8 = 777 at neutron temperature of disk. 
A235 = atomic weight of U-235 

A238 = atomic weight of U-238 

AC5 = oat neutron baccenties of 340°K 

AC8 = Tt at neutron temperature of 340° K 

AN = Avogadro's number 

apa = (— : 

AU02 = atomic weight of uranium dioxide 

B= buckling 

CCD = C for a fuel disk with change in moderator mass 


CD = C for a fuel disk 


DD = D for a fuel disk 





DEAP = B24 
pELG = 30. 
DELD = 3D~ 


DELEA = Ske 
DELEF = 2, 


DELH disk thickness 


3P : 


DELP 





DELPLC 


change in density of polyethylene within the fuel core 


DELPLD = change in density of polyethylene within a fuel disk 


DELSAS = ‘%0,*** 
DELSA8 = $0" 
DELSF = 24) 
DELTA = A&A 


DENP = density of polyethylene in a disk 

DENPLD = density of polyethylene in a disk with change in moderator mass 
DENUO2 = density of uranium dioxide in a disk 

DENUS = density of U-235 in a disk 


DENUS 


density of U-238 in a disk 


DDD = D for a fuel disk with change in moderator mass 


EAP = ed 

EAU = be 

EO= — 8 
BPA = de 


ERI = effective resonance integral 
© #e ? 
ESP = dg 
EST= By 
reas v 
_Esu0z2 = &s 
GA5 = non-1/v absorption factor for U-235 


GAB 


non-1/v absorption factor for U-238 
GF = non-i/w fission factor for U-235 


H = distance of nearest edge of fuel disk above or below the zero 
ordinate of the Z-axis of the cylindrical fuel core 


HH = height of fuel core 
HP = extrapolated height of fuel core 
HR = equivalent height of reflected fuel core 








PD = P for a disk 
PM = mass of polyethylene in a fuel disk 
PPD = P for a disk with change in moderator mass . 
RHOD = reactivity of fuel disk 
RHOR1 = reactivity of fuel core for bare core geometry model 
RHOR2 = reactivity of fuel core for extrapolated dimensions model 
RHOR3 = reactivity of fuel core for reflected core model 
_ SA50 = Go at neutron temperature of 293° K 
SA80 = oq at neutron temperature of 293° K 


= 265 
SABARS = Wa 


= 238 
SABARS8 = 


SFBAR = 4 

SFO = {; at neutron temperature of 293° K- 
SIBAR = & 

sip= §° 

stuo2 = §&° 

SN = =, 

spa = (-)! 

suoz = ,° 

TEFF = effective neutron temperature 


TM = moderator temperature 


TO= @, 


VD = volume of disk 
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END 


fe BCORE (DELP »DELC,NELEF, DEL EA, DEL Dy Ny RELELC, Hy DELH,P@ HOR 1 

aoe H(10) ,DELH(10) »RHON(10) 

8: «142 ADELH (ND /HH} 40. 5*SIN(6.294% (HON) 4DELH(N )) /HHI-0 CEXSIN( 6.28 
42 aDELHIND /HH)~0. 5*SIN(6 2.284% (HIND 4D ELHIN ))/SHHD 40. 5*SIN( 6228 


/3,. MABE eae ees 34O*XDEL C415. 37* DEL EF-8. B44eDELFA~ (0. 
9*Y/Q) & DO) 


oe et 
wwe 5 


10 


o “O0ON” = 
— 
Suse 
Nees 


4X, 8 REACTIVITY =',F1C.4) 


2H MO we te 
e 
~ 


Cxyuo Omme ae 


oo 
fasta 
+0 


+#RHON (I) 


mxXie 


i 
Oo 
A ed 
MVODWDIDAaANswWIH< 


ZNGQTOTDANO DH =z # il 
woe Ve DIT mye we 


CADWREOZEINGAwW 


94429 O4e0enNOTs 
Zz it 


~~ 


JBROUTINE EXTDIM (DELP,DELC,DELEF, DEL FA, DEL DN, DELPLC yy HyDFLHyS HOF 
ENSTON DELH(10) gHO10) »RHON(10) 


~ 


P)+0.25¥*SIN(6.284*HH/HP) 
Ghee tLe 264% (H(N)+40EL HUN) )/HP )-9.5*S ING 6, 24x 


ptt 
“Nee ols tar CONM 


NDVSLl= i ~ i WoL 


on Kee 


M 
36. 
3 
N 
3. 
N 
C 


~ 
Re til e~— en 
ir 


bat SD pet OS oe, 


~ 
© 


91Xo1254X_, ' PEACTIVITY =',F10.6) 


ZeqANOwe we 6 Ow! 
e Perrin ew erie ° 
Mame TSS Bou 
me er ORONO 


92 


1¢ 


20 


30 


a) 


ph 


RETURN 


Oo RHOR2=C.9 


D0 30 J=1510_ 
RHCRO= RHOK2+RHAD(S) | 


SUBROUTINE RFFLCT (DELP,DELC,DELEF, DEL FA, DELO yN ep DELPLCy Hy DELH,RHOAR 
DIMENSTON DELHILO) HELO) RHOD(10) 


HH=9.5 

HR=12,25 

patna DELHCN) JHR4+O5¥S ING. 284% (HIN) FDEL HEN) I/HR I-0. SRSING 6 0 234 
Se SESIN(GE. 284% (HIN) +DELHIN) D/HR +9. S*SING GS. 234 
G=3.142*HH/ (2. cot ttcene ners 28 4*UHH/HR ) 

RHOD (NDE CF /(2 245.0) «(1 205 NELB4+1. 349%XDFLC+15. ST*DELEF- 2. RLGENELEL= 
WRITE (6910) NyRHODIN) | 

EQRMAT(B5X,_, !' DISK, IX,T2,4X_  ' REACTIVITY =", F16.45) 

TF(NJE9210) GO TO 20 

RETURN - 

RHOR3=0.09 

99 3C Jj=1,10 

RHCR3=RHOR3B4+RHOD (I) 

CONTINUE 

RE TURN S 

END 


APPENDIX F_ 
COMPUTER PROGRAM 2 


Ao COMPUTER. NOTATION. 

The majority of the symbols appearing in this computer program 
are the same as in Computer Program 1 and have the same meaning. 
Therefore, reference shouta be made to Appendix E with the exception. 


of the following terms... 


density of polyethylene in new fuel disk 


DENPND = 
DNUSND = density of U-235 in new fuel disk 
DNU8SND = density of U-238 in new fuel disk 
DUO2ND = density of uranium dioxide in new fue? disk 
_ RHOMAX © = value of excess. reactivity of the core etter interchanging 


old and new disk 
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APPENDIX G 
RESULTS OF COMPUTER PROGRAM 1 


TABLE 6 i 


Change in Moderator Core Mass = 10. gm 





Reactivity Reactivity 
Reactivity _ (Extrapolated (Reflected 
Disk (Bare Core) Dimensions) Core) 
1 ~-0.000000 | | -0.000000 -0.000000 | 
2 -0.000001 -0.000001 -0.000003 | 
3 -0.000002 -0.000002 -0.000004 
4 -0.000006 -0.000007 -0.000012 : 
5 ~0.000028 -0.000034 -0.000045 f 
«6 -0.000050 -0.000059 -0.000068 , 
7 -0.000112 -0.000129 - 0.000137 ‘ 
8 -0.000116 -0.000134 -0.000140 | 
9 ~-0.000079 -0.000091 ~0.000097 ! 
10 -0.000016 -0.000019 | 


~0.000025 


Reactivity of core (bare core solution) = -0.000410 


Reactivity of core (extrapolated dimensions solution) = ~0.000476 


Reactivity of core (reflected core solution) = -0.000530 
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TABLE 7 


« Change in Moderator Core Mass = 20. gm 
| | scaeeicies: Reactivity 
a Reactivity (Extrapolated (Reflected 
Disk ‘(Bare Core) . Dimensions) Core) 
1 -0.000000 -0 000000 ~ 0.000001 
2 -0.000002 -0 000003 -0.000006 
3 ~Q.000004 ~0,000005 -0.000010 
4 -0.000015 -0.000018 -0.000029 
5 -0,000051 -0,000061 -0.000081 
6 -0.000093 -0.000110 -0.000127 
7 -0.000229 -0.000265. -0.000280 
8 ~0,000236 = ,.000272 -0.000284 
. 9 -0.000154 -0.000178 -0.000191 
10 -0.000037 


~0.000032 


~0.000048 


Reactivity of core (bare core solution) = -0.000818 


Reactivity of core (extrapolated dimensions solution) = =0,000949 


Reactivity of core (reflected core solution) = ~0 001058 
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‘TABLE 8 


Change in Moderator Core Mass = 30. gm 





Reactivity Reactivity 
Reactivity (Ext rapolated (Reflected 
Disk (Bare Core) Dimensions) Core) 

1 -0.000000 -0.000000 -0.000001 
2 -0.000004 ~.0000004 -0.000009 
3 ~0.000007 ~0.000008 -0.000016 
4 ~0.000024 -0.000028 -0.000047 
5 ~0.000075 -0.000088 ~0.000118 
6 ~0.000137 -0.000161 -0.000187 
7 -0.000347 -0.000401 -0.000424 
8 -0.000356. -0.000411 -0.000428 
9 -0.000230 -0 000266 -=0..000285 
10 ~0.000047 -0.000055 -0.000071 


Reactivity of core (bare core solution) = -0.001227 


Reactivity of core (extrapolated dimensions solution) = -0.001423 


Reactivity of core (reflected core solution) = -0,001587 
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TABLE 9 


Change in Moderator Core Mass = 40. gm 


. Reactivity. Reactivity 
Reactivity (Extrapolated (Reflected 
Disk (Bare Core) Dimensions) Core) 
1 -0.000000 0.000000 -0.000002 
2 0.000005 0.000005 0.000012 
3 -0.000010 -0,000011 -0.000022 
4 -0.000033 0.000039 ~0.000065 
5 -0.000098 | -0.000116 -0.000155 
6 -0.000181 ~0.000213 0.000246 
7 0.000465 -0.000537 -0.000568 
. 8 -0.000476 -0.000550 -0.000574 
9 -0.000306 0.000354 -0.000379 
. 10 0.000073 -0.000095 


-0.000062 


Reactivity of core (bare core solution) = -0.001638 


Reactivity of core (extrapolated dimensions solution) = -0,.001899 


Reactivity of core (reflected core solution) = -0.002119 
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TABLE 10 


Change in Moderator Core Mass = 50. gm 





Reactivity Reactivity a 

Reactivity (Extrapolated (Reflected 

Disk (Bare Core) Dimensions) Core) 
1 -0.000001 -0.000001 '-0.000002 
2 -0.000006 -0.000007 -0.000016 
3 -0.000013 ~0.000014 -0.000029 
4 -0.000042 -0.000050 =0.000084 
5 -0.000121 -0.000143 -0.000192 
6 ~0,000225 -0..000265 -0..000306 
7 -0.000583 -0.000674 -0.000713 
8 ~0.000597 -0.000690 0.000719 
“ug 0.000382 -0.000442 -0.000472 
10 ~0.000077 -0.000090 -0.000118 


Reactivity of core (bare core solution) = -0,002049 


Reactivity of core (extrapolated dimensions solution) = -0.002376 


Reactivity of core (reflected core solution) = -0.002651 


86 








APPENDIX H 
RESULTS OF COMPUTER PROGRAM 2 


TABLE 11 


Reactivity Change Resulting From 
Interchanging New Disk With Disk 4 


Maximum 
; _ Excess — 
Model Reactivity Reactivity ~ 
Bare Core 0.001582 0.004392. 
Extrapolated Dimensions 0.002006 = 0.004816 
Reflected Core 3 0.002957 — -0,005767 
TABLE 12 
Reactivity Change Resulting From 
Interchanging New Disk With Disk 5 | 
Maximum 
. "os Excess . 
Model . Reactivity Reactivity — . 
Bare Core 0.002490 0.005300 
Extrapolated . 0.003128 0.005938 
Reflected Core -0.003833 0.006643 
TABLE 13 
Reactivity Change Resulting From 
interchanging New Disk With Disk 6 
Maximum 
Excess 
Model Reactivity Reactivity 
Bare Core 0.003671 0.006481 
Extrapolated Dimensions 0.004541 6.007351 
Reflected Core 0.004897 te 0.007707 
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